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Abstract 

Cycloaddition  reactions  between  1,3-butadiene  and  the  C-terminated  SiC(l  0  0)-c(2x2)  surface  have  been  addressed 
using  quantum-chemical  methods.  The  c(2x2)  structure  consists  of  — C=C—  bridges  between  underlayer  Si  atoms 
which  themselves  form  Si— Si  bonds.  Of  various  possible  reaction  products,  the  one  formed  by  a  [2  +  4]  reaction  with  the 
— C=C—  bridge  (giving  a  species  resembling  1,4-cyclohexadiene)  is  the  lowest  in  energy.  Density  functional  calculations 
for  the  bare  c(2  x  2)  surface,  using  a  cluster  model  with  mechanical  embedding,  gave  good  agreement  with  structural 
parameters  obtained  in  previous  fully  ab  initio  studies.  Similar  calculations  for  the  cycloaddition  product  and  for  the 
transition  state  gave  a  reaction  energy  of -50.3  kcal/mol  and  an  activation  energy  of  +6.1  kcal/mol  to  form  a  planar  ring 
structure  lying  normal  to  the  surface.  Detailed  results  for  the  frequency  and  infrared  polarization  behavior  of  adsorbate 
vibrational  modes  have  also  been  obtained. 

©  2003  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Cycloaddition  reactions  of  alkenes  and  dienes 
with  weakly  Ti-bonded  dimers  on  the  (1  00)-(2x  1) 
surfaces  of  the  column-IV  semiconductors  C  (di¬ 
amond),  Si  and  Ge  have  received  much  attention 
[1-3].  Such  reactions  provide  an  important  class  of 
methods  for  fabricating  “functionalized”  semi¬ 
conductor  surfaces  with  potential  applications  in 
chemical  sensors  and  in  hybrid  electronic  devices. 
Ab  initio  quantum-chemical  modeling  [4-12]  has 
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been  a  powerful  tool  in  understanding  structure 
and  reactivity  in  these  cycloaddition  processes. 

In  this  work,  we  consider  cycloaddition  reac¬ 
tions  between  1,3-butadiene  (the  prototypical  di¬ 
ene)  and  the  (100)  surface  of  the  IV-IV  compound 
semiconductor  cubic  SiC.  The  SiC(lOO)  surface 
exhibits  several  reconstructions,  depending  on 
stoichiometry,  which  have  been  reviewed  recently 
from  both  experimental  [13,14]  and  theoretical  [15] 
perspectives.  Here  we  focus  on  the  very  unusual  C- 
terminated  c(2  x  2)  structure  consisting  (see  below) 
of  essentially  triple-bonded  —  C=C—  bridges  be¬ 
tween  Si  atoms  which  themselves  form  Si— Si  cj- 
bonds.  The  resulting  surface  has  no  dangling  bonds 
and  no  intrinsic  surface  states  in  the  band  gap. 
Previous  theoretical  studies  [16,17]  have  addressed 
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the  chemisorption  of  small  molecules  on  SiC(  1  0  0) 
surfaces;  although,  to  our  knowledge,  the  particu¬ 
lar  behavior  of  the  c(2  x  2)  in  this  regard  has  been 
largely  unexplored. 


2.  Results  and  discussion 

2.1.  Possible  adsorbate  structures 

Fig.  1  shows  a  schematic  diagram  of  the 
SiC(l  00)-c(2x2)  surface.  Fig.  2  illustrates  six 
different  adsorbate  structures  which  could,  hypo¬ 
thetically,  be  formed  by  reaction  with  1,3-butadi¬ 
ene.  The  bare  c(2x2)  surface  will  be  examined  in 
detail  below. 

The  first  step  was  to  evaluate  the  relative  ener¬ 
gies  of  the  various  adsorbate  structures.  This 
was  done  by  constructing  the  cluster  model  for 
the  bare  c(2  x  2)  shown  in  Fig.  3a  and  using  either 
a  semiempirical  or  a  molecular  mechanics  (MM) 
approach  to  optimize  the  geometries  of  the 


Fig.  1.  Schematic  model  for  the  SiC(l  00)-c(2x2)  surface 
[13,15]  viewed  along  the  surface  normal  (a)  and  in  the  surface 
plane  (b).  The  cross-hatched  circles  are  the  — C=C—  bridges 
lying  above  the  plane.  The  open  circles  the  first-underlayer  Si 
atoms  lying  in  the  plane,  and  the  filled  circles  the  second- 
underlayer  C  atoms  lying  below  the  plane.  Dashed  lines  show 
the  bonds  between  Si  atoms,  and  the  dashed  square  shows  the 
c(2x2)  surface  unit  cell. 


Fig.  2.  Schematic  models  for  the  six  possible  1,3-butadiene  re¬ 
action  products  considered  in  the  present  work.  The  cis  and 
trans  molecular  structures  are  also  shown.  The  small  filled  circles 
in  (c)  and  (d)  indicate  unpaired  electrons.  Heavy  lines  indicate 
>C=C<  tt-bonds,  and  extra-heavy  lines  show  — C=C—  7>bonds. 
Dashed  lines  show  Si— Si  a-bonds.  For  clarity  the  H  atoms  and 
most  of  the  C  atoms  are  not  shown  explicitly.  The  structures  are 
idealized  models  and  have  not  been  optimized. 


Fig.  3.  Schematic  diagram  of  the  cluster  model  for  the  c(2x2) 
surface,  (a)  shows  the  entire  SfeCjsH^  cluster,  and  (b)  shows 
the  "high-level”  part  treated  via  DFT.  The  remainder  of  (a) 
forms  the  “low-level”  part  treated  using  RHF/PM3  or  UFF 
methods.  The  various  atom  types,  including  the  S  link  atoms, 
are  indicated.  The  structure  shown  is  the  result  of  the  geometry 
optimization. 
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cluster  +  adsorbate  (bonded  as  in  Fig.  2).  For  the 
former,  restricted  Flartree  Fock  (RHF)  calcula¬ 
tions  were  done  at  the  PM3  level  [18].  Badziag  [19] 
has  performed  a  purely  semiempirical  treatment  of 
the  bare  c(2x2)  surface  at  the  RHF/PM3  level, 
using  a  periodic  cluster  model,  and  obtained  good 
agreement  with  ab  initio  geometries  (see  below). 
For  MM,  the  Universal  Force  Field  (UFF)  [20]  in 
the  Gaussian98  program  suite  [21]  was  used. 

Structures  based  on  reaction  across  the  Si— C 
bond  (Fig.  2c,  d)  are  unfavorable  because  they 
result  in  unpairing  of  electrons  in  the  — C=C—  and 
Si— Si  bonds  of  the  bare  surface.  These  were  not 
further  considered  in  the  present  calculations.  Of 
the  remaining  four  structures,  both  methods 
(RHF/PM3  and  UFF)  agree  that  the  product 
formed  by  [2  +  4]  cycloaddition  at  the  — C=C— 
bridge  (Fig.  2a)  is  the  most  stable.  The  RHF/PM3 
result  places  this  at  ~38  kcal/mol  below  the  sec¬ 
ond-most  stable  structure,  the  [2  +  2]  product  (Fig. 
2b).  A  more  accurate  value  for  this  energy  differ¬ 
ence  is  given  below.  The  prediction  of  the  [2  +  4] 
(vs.  [2  +  2])  cycloaddition  product  as  being  the 
more  stable  is  consistent  with  ab  initio  results  for 
the  reaction  of  1,3-butadiene  with  the  (1  00)-(2x  1) 
surfaces  1  of  C  [4],  Si  [5,8]  and  Ge  [8],  It  is  also 
noted  that  the  [2  +  2]  reaction  is  symmetry-for¬ 
bidden,  according  to  the  Woodward-Hoffman 
rules,  for  a  symmetric  — C=C—  bridge  as  in  Fig. 
lb.  However,  on  the  real  surface  (see  below)  the 
bridge  is  asymmetric,  which  makes  the  [2  +  2]  for¬ 
mally  allowed  (although  still  energetically  unfa¬ 
vorable).  In  the  following,  attention  will  focus 
mainly  on  the  [2  +  4]  reaction,  with  some  further 
mention  of  the  [2  +  2]. 

2.2.  The  bare  c(2  x2)  surface 

A  cluster  model  was  constructed  for  the  bare 
c(2x2)  surface  as  a  basis  for  treating  the  adsor¬ 
bate.  This  was  done  using  a  mechanical-embed¬ 
ding  (ME)  approach  [22],  as  implemented  in  the 
Gaussian98  program  suite  [21],  which  has  been 
applied  successfully  to  surface  problems  [23],  The 


1  The  Ge  results  in  [8]  are  for  Ge-Ge  dimers  on  a  Si(  1  0  0) 
surface. 


ME  technique  used  was  that  recently  developed  by 
Kruger  and  Sax  [24]  involving  divalent  atoms  (Be, 
O  or  S),  rather  than  H,  as  link  atoms.  (Link  atoms 
replace  cluster  atoms  during  the  high-level  treat¬ 
ment  of  the  chemically-significant  sub-section  of 
the  cluster.)  Fig.  3a  shows  the  entire  cluster,  and 
Fig.  3b  shows  the  part  that  was  treated  at  the  high 
level,  together  with  the  S  link  atoms.  The  surface 
Si— Si  bonds  are  an  integral  part  of  the  c(2x2) 
structure,  and  this  sub-cluster  is  the  smallest  that 
can  possibly  be  expected  to  include  these  ade¬ 
quately.  Forming  the  high-level  sub-section  re¬ 
quires  breaking  two  Si— C  bonds.  Initial  attempts 
to  use  H  as  the  link  atom,  which  then  gives  two 
Si— H  bonds  per  high-level  Si  atom,  led  to  diffi¬ 
culties  in  geometry  optimization  which  were 
eliminated  by  the  present  approach. 

The  low-level  part  of  the  ME  calculation  was 
done  using  either  RHF/PM3  or  UFF  methods. 
The  UFF  atom  types  [20]  were  declared  to  be 
“H_”  and  “Si3”  for  H  and  Si,  “C_3”  and  “C_l”, 
respectively,  for  sub-surface  C  and  — C=C—  and 
“S_2”  for  the  S  link  atoms.  The  high-level  treat¬ 
ment  of  the  sub-cluster  (Fig.  3b)  was  done  using 
density  functional  theory  (DFT)  with  the  B3LYP 
functional,  and  no  constraints  were  applied  in  the 
geometry  optimization.  In  the  following,  calcula¬ 
tions  will  be  termed  “(A/B:C)”  where  “A”  is  the 
high-level  method  (usually  B3LYP),  “B”  the  high- 
level  basis  set  (usually  6-3 1G*)  and  “C”  the  low- 
level  method  (either  PM3,  meaning  RHF/PM3,  or 
UFF). 

Table  1  (cf.  Fig.  4a)  compares  structural  results 
with  those  obtained  using  fully  ab  initio  slab- 
model  treatments  [25-28]  and  with  the  available 
experimental  data  [29-31].  In  the  following,  only 
the  central  — C=C—  bridge  will  be  considered,  in 
order  to  avoid  ambiguities  due  to  edge  effects.  All 
calculations  give  a  — C=C—  distance  of  about  1.21 
A,  close  to  the  corresponding  value  [32]  in  the 
molecule  H3SiC=GSiH3  (1.216  A).  The  Si-Si 
distance  found  in  previous,  fully  ab  initio  work 
[25-28]  is  close  to  that  for  dimers  on  the  Si(l  0  0)- 
(2xl)H  monohydride  surface  (2.39  A  [5]).  In 
either  system,  there  are  no  unpaired  dangling 
orbitals  on  surface  Si  atoms.  The  (B3LYP/ 
6-31G*:PM3)  value  obtained  here  (f  =  2.42  A)  is 
in  good  agreement  with  these  results.  The  (B3LYP/ 
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Table  1 


Structural  results  for  the  ME  cluster  model  (“this  work”)  vs.  fully  ab  initio  treatments  of  the  SiC(l  00)-c(2x2)  surface 


Distance 

Yan  [25] 

Kackell  [26] 

Sabisch  [27] 

Catellani  [28] 

This  work3 

Expt. 

d\ 

1.22 

1.23 

1.22 

1.23 

1.23  (1.23) 

1.22b 

1.31° 

d2 

1.87 

1.82 

1.83 

1.81 

1.85  (1.84) 

1.84b 

1.93° 

dz 

2.38 

2.38 

2.40 

2.37 

2.42  (2.53) 

2.70b 

2.71° 

hz\ 

1.30 

1.32 

1.32 

1.32  (1.30) 

1.60c 

8z2 

1.02 

1.04 

1.08  (1.12) 

1.07° 

All  distances  are  in  Angstroms.  The  terms  are  defined  in  Fig.  4a. 

“Values  given  are  the  (B3LYP/6-31G*:PM3)  results.  Those  in  parentheses  are  from  the  (B3LYP/6-31G*:UFF)  calculation. 
b  Photoelectron  diffraction  data  [29,30]  for  a  sample  prepared  by  thermal  desorption  of  Si  from  a  Si-terminated  surface. 

“Low  energy  electron  diffraction  data  [31]  for  a  sample  prepared  by  desorption  of  Si.  Surfaces  prepared  by  exposure  to  C2H4  at  high 
temperature  gave  somewhat  different  distances,  possibly  due  to  the  presence  of  adsorbed  H. 


6-31G*:UFF)  result,  while  somewhat  greater 
(di  =  2.53  A),  is  still  much  less  than  the  Si— Si 
distance  on  the  ideally-terminated  SiC(lOO)  sur¬ 
face,  3.08  A.  Thus,  the  ME  treatment  of  the  cluster 
shown  in  Fig.  3  appears  to  be  comparable  to  fully 
ab  initio  models  in  describing  — C=C—  bridges  on 
the  c(2x2)  surface.  None  of  the  results,  including 
those  obtained  here,  show  the  experimentally- 
observed  tilting  (“buckling”)  of  the  — C=C— 
bridge  [33,34].  Recent  work  [35]  indicates  that  the 
tilting  may  be  due  to  stress  along  the  — C=C— 
bridge,  which  was  not  considered  in  the  models 
included  in  Table  1.  Also,  the  calculated  Si-Si 
distance  (c/3)  is,  in  all  cases,  significantly  shorter 
than  the  experimental  value. 

Unrestricted  (UB3LYP)  calculations  were  done 
to  test  for  the  possibility  that  the  Si— Si  bond 
might  be  so  weak  that  unpaired  spin  density  would 
appear  on  surface  Si  atoms.  However,  the  calcu¬ 
lation  always  converged  to  a  closed-shell  singlet. 
The  vibrational  frequencies  of  the  optimized 
structure  were  checked  for  the  absence  of  imagi¬ 
nary  values  which  would  indicate  a  geometric  in¬ 
stability.  Of  the  possible  divalent  link  atoms  Be,  O 
and  S  [24],  preliminary  (UHF/3-21G:UFF)  calcu¬ 
lations  found  that  only  S  gave  good  convergence 
in  geometry  optimization.  The  optimized  geometry 
was  remarkably  insensitive  to  the  choice  of  high- 
level  method  and  basis  set.  A  comparison  of 
(UHF/3-21G:UFF)  and  (B3LYP/6-31G*:UFF) 
results  showed  differences  in  d\ ,  and  of  2% 
and  in  8z,  (§r2)  of  3%  (7%).  Table  1  also  shows 


that,  except  as  noted  above,  the  structure  is  es¬ 
sentially  independent  of  the  low-level  method. 

2.3.  The  [2  +  4]  cycloaddition  product:  structure 
and  reaction  energy 

The  cluster  model  was  augmented  to  include  the 
[2  +  4]  product  at  the  central  — C=C—  bridge  (Fig. 
4b),  and  Table  2  gives  the  optimized  structural 
parameters.  The  adsorbate  resembles  1,4-cyclo- 
hexadiene  (Fig.  4c)  with  the  plane  of  the  six- 
member  ring  lying  normal  to  the  surface.  Both 
theoretical  [36]  and  spectroscopic  [37]  results  in¬ 
dicate  that  the  free  1 ,4-cyclohexadiene  ring  is  pla¬ 
nar,  and  the  same  was  found  here  for  the 
cycloaddition  product.  The  absence  of  imaginary 
vibrational  frequencies  indicates  that  this  is  a  sta¬ 
ble  structure.  The  d\  and  d(,  distances  are  typical  of 
a  C=C  71-bond,  and  the  d4  and  d4  values  are  typical 
of  a  C— C  a-bond.  The  C— H  distances  are  1.088  A 
for  C— C — H  and  1.103  A  for  CH2.  These  value 
and  ds  and  d6  are  all  close  (within  ~0.01  A)  to  the 
corresponding  values  reported  in  a  recent  theo¬ 
retical  study  of  molecular  1 ,4-cyclohexadiene  [36]. 
Distances  involving  only  substrate  atoms  (i.e.,  d4 
and  8z2)  are  essentially  unaffected  by  adsorption 
(compare  corresponding  methods  in  Tables  1  and 
2).  This  suggests  that  underlayer  atoms,  particu¬ 
larly  the  Si— Si  bond,  are  not  strongly  involved  in 
chemisorption  at  the  — C=C—  site  and  supports 
the  use  of  the  small  sub-cluster  (Fig.  3b). 
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Fig.  4.  (a)  Schematic  diagram  showing  the  distances  (cf.  Table 
1)  defining  the  c(2x2)  surface  structure.  The  extra-heavy  C— C 
lines  indicate  the  — C=C —  7t-bonds.  (b)  Schematic  diagram  of 
the  [2  +  4]  cycloaddition  reaction  product  for  1,3-butadiene  and 
the  c(2x2)  surface  (cf.  Table  2).  The  heavy  lines  indicate 
>C=C<  tt-bonds.  For  clarity,  the  ring  C  atoms  and  the  two  FI 
atoms  on  one  of  the  C  atoms  are  omitted,  (c)  The  free  1,4-cy- 
clohexadiene  molecule.  The  different  coordinate  systems  for  the 
symmetry  analysis  of  (b)  and  (c)  are  indicated. 

To  obtain  a  value  for  the  energy  of  reaction 
(A Er),  defined  as  the  energy  difference  between  the 
product  and  the  reagents  at  0  K,  similar  calcula¬ 
tions  were  done  for  free  1,3-butadiene.  Table  3 
compares  the  result,  A ER  =  -54.8  eV,  with  those 
for  the  same  reaction  on  other  surfaces.  For  the 
present  results,  the  electronic  energies  were  cor¬ 
rected  for  the  zero-point  (ZP)  vibrational  energies, 
and  the  cis  form  of  1,3-butadiene  was  used  as  the 
reagent  (as  required  by  the  geometry  of  the  [2  +  4] 
product).  The  trans  form  was  computed  to  be 
lower  in  energy  by  3.5  kcal/mol;  hence,  AEr  would 


Table  2 

Structural  parameters  (A),  defined  in  Fig.  4b,  for  the  reaction 
product  and  for  the  transition  state 


Distance 

Product 

Transition  state 

d\ 

1.36  (1.36) 

(1.26) 

d2 

1.92  (1.91) 

(1.84) 

d2 

2.49  (2.55) 

(2.54) 

d4 

1.53  (1.53) 

(2.43) 

d5 

1.50  (1.50) 

(1.37) 

de 

1.33  (1.33) 

(1.43) 

hz\ 

1.57  (1.55) 

(1.35) 

8z2 

1.06  (1.11) 

(1.12) 

Values  in  parentheses  are  from  (B3LYP/6-31G*:UFF)  calcu¬ 
lations.  The  others  are  (B3LYP/6-31G*:PM3)  results. 


be  less  exothermic  by  this  amount  if  the  trans  form 
were  used  as  the  reagent.  It  is  likely  that  strain 
makes  a  significant  contribution  to  AFr.  For  the 
bare  c(2  x  2)  surface  the  computed  Si— C=C  angle 
is  135°,  far  removed  from  the  ideal  sp-hybridized 
value  of  1 80°;  whereas,  in  the  reaction  product  the 
Si— C=C  angle  is  126°,  close  to  the  ideal  sp2  value 
of  120°. 

Some  of  the  A ER  values  in  Table  3  are  for  dif¬ 
ferent  basis  sets,  but  the  differences  for  the  various 
semiconductors  are  greater  than  the  expected  basis- 
set  dependencies.  For  example,  A FR  for  the  [2  +  4] 
cycloaddition  reaction  of  1,3-cyclohexadiene  with 
Si(l  0  0)-(2x  1)  [5]  is  <3%  more  exothermic  for  6- 
311G**  than  for  6-3 1G*  basis  sets.  In  the  present 
case,  the  basis-set  dependence  was  examined  using 
single-point  (B3LYP/6-311+G*:PM3)  calcula¬ 
tions  for  structures  optimized  at  the  (B3LYP/ 
6-31G*:PM3)  level.  The  ZP  corrections  at  the 
(B3LYP/6-31G*:PM3)  level  were  also  used.  The 
result  is  AFR  =  -50.3  kcal/mol,  ~8%  less  exo¬ 
thermic  than  the  all-(B3LYP/6-31G*:PM3)  result. 
It  is  speculated  that  the  larger  basis-set  dependence 
in  this  case  may  indicate  that  including  a  diffuse 
function  improves  the  description  of  the  weak  (i.e., 
long)  Si-Si  bond.  With  B3LYP/6-31G*  as  the 
high-level  method,  little  or  no  difference  is  noted 
between  results  obtained  with  RHF/PM3  or  UFF 
as  the  low-level  method.  This  will  be  useful  in 
treating  the  transition  state  (see  below). 

Similar  calculations  were  also  performed  for  the 
[2  +  2]  reaction  product.  The  ZP-corrected 
(B3LYP/6-31G*:PM3)  result  places  this  species 
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Table  3 

Comparison  of  reaction  (A f?R)  and  activation  (AEa)  energies  computed  for  the  [2  +  4]  cycloaddition  reaction  of  1,3-butadiene  with 
various  surfaces 


Substrate 

aer 

AEa 

C(1 0  0)-(2  x  1 ) 

-76.2“ 

+5.2 

-84.3  to  -85. 5b 

+2.7  to  +1.9 

Si(l  0  0)-(2  x  1) 

-67. 5C 

—6 1 .8d 

Ge(l  00)-(2x  l)e 

-45. 8d 

SiC(100)-c(2x2) 

-54.8  (-56.4/ 

(+4.7) 

-50.38 

(+6.1) 

All  values  are  in  kcal/mol,  and  negative  values  indicate  an  exothermic  process.  All  studies  used  DFT  with  the  B3LYP  functional.  The 
footnotes  give  the  basis  sets,  the  form  (“ cis ”  or  “trails")  of  1,3-butadiene  used  as  the  reagent  (see  text),  and  whether  it  was  explicitly 
stated  that  electronic  energies  were  corrected  for  zero-point  (ZP)  vibrational  energy. 

“Ref.  [4]  (6-31+G*;  cis;  ZP-corrected). 

bRef.  [11]  (6-31G*;  cis;  ZP-corrected).  Results  depend  slightly  on  cluster  size. 

°  Ref.  [5]  (6-3 1G*;  cis). 

dRef.  [8]  (6-31G*  for  terminating  H  and  bulk  Si,  6-31 1+G**  for  surface  Si  or  Ge  and  for  1,3-butadiene;  trans). 

“These  results  are  for  Ge-Ge  dimers  on  a  Si(l  00)  surface. 

fPresent  work  [(B3LYP/6-31G*:PM3);  cis;  ZP-corrected],  Values  in  parentheses  are  (B3LYP/6-31G+UFF)  results.  The  trans  form 
of  1,3-butadiene  was  computed  to  he  3.5  kcal/mol  (ZP-corrected)  lower  than  the  cis  at  the  B3LYP/6-31G*  level. 

SA£R  as  in  f  but  for  (B3LYP/6-311+G*:PM3)  single-point  calculations  using  the  (B3LYP/6-31G*:PM3)  structures  and  ZP-cor- 
rections  obtained  in  f.  The  A EA  value  (in  parentheses)  is  the  (B3LYP/6-31 1+G*:UFF)  single-point  result  obtained  using  the  (B3LYP/6- 
31G*:UFF)  transition  state  structures  and  ZP-corrections  in  f. 

23.5  kcal/mol  above  the  [2  +  4]  product.  The  pre¬ 
sent  result  is  close  to  the  [2  +  2]-[2  +  4]  difference 
computed  for  1,3-butadiene  reacting  with  Si  and 
diamond  (10  0)-(2xl).  For  Si,  values  (not  cis- 
trans  corrected)  of  25.6  [5]  and  24.8  kcal/mol  [8] 
have  been  found.  For  diamond,  a  value  of  28.2 
kcal/mol  [4]  has  been  reported.  For  both  Si(10  0) 

[38]  and  C(10  0)  [4],  analysis  of  infrared  (IR) 
spectroscopic  data  indicates  that  the  reaction  is 
predominantly  [2  +  4],  and  the  same  is  expected  for 
SiC(l  00)-c(2x2). 

2.4.  Vibrational  spectra 

The  frequencies  of  adsorbate  vibrational  modes 
are  compared  in  Table  4  with  corresponding  val¬ 
ues  from  IR  [39,40]  and  Raman  [37]  data  for  liquid 
1,4-cyclohexadiene.  The  calculated  modes  are  all 
localized  within  the  adsorbate;  although,  those 
below  ~1100  cm-1  exhibit  small  contributions 
from  other  atoms.  Hence,  the  frequencies  are  ac¬ 
curate  at  the  B3LYP/6-31G*  level  [22]  and  have 
been  scaled  accordingly  [41].  Some  free-molecule 
modes  will  of  course  be  strongly  perturbed  in  the 
adsorbate.  Results  for  RHF/PM3  or  UFF  (not 


shown)  as  the  low-level  method  differed  by  ^4 
cm-1  except  for  the  stretching  mode  of  the  C=C 
bonded  to  Si  (1561  cm-1  in  Table  4)  which  was  25 
cnr1  lower  for  UFF.  Vibrational  modes  were  also 
obtained  for  the  [2  +  2]  product  but  are  not  given 
here. 

The  relevant  surface  experiment  [1]  is  typically 
done  using  IR  spectroscopy  in  the  multiple  internal 
reflection  (MIR)  configuration  2  which  provides 
polarization  data.  In  principle,  high-resolution 
electron  energy  loss  spectroscopy  could  also  be 
used.  However,  for  a  partially  ionic  material  such 
as  SiC,  detection  of  weak  adsorbate  modes  is 
complicated  by  the  presence  of  strong  Fuchs- 
Kliewer  phonon  losses  [42], 

Issues  related  to  symmetry  are  now  considered. 
The  point-group  of  the  free,  planar  1 ,4-cyclohexa- 
diene  molecule  is  D2h,  and  that  of  the  adsorbate 
(Fig.  4b)  is  C2v  with  ( x,z )  and  ( y,z )  as  mirror 


2  Only  adsorbate  modes  calculated  to  lie  above  1000  cm-1 
are  given  in  Table  4.  For  an  MIR  experiment  using  as  a  sample 
a  thin  film  of  SiC  grown  on  a  Si  substrate,  this  would  be  the 
approximate  lower  limit  of  the  accessible  range  imposed  by  the 
transmission  of  the  Si  MIR  element. 
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Table  4 

Computed  vibrational  frequencies  (cm-1)  for  the  cycloaddition  product  (“adsorbate”)  together  with  corresponding  IR  and  Raman 
values  for  liquid  1 ,4-cyclohexadiene 


Modea 

Adsorbate15 

IR  [39]c 

Raman  [37] 

=CH  sym.  str.d 

3061  (z;  100) 

3031  (Ag,  B2g)e 

=CH  asym.  str.d 

3040  (y;  12) 

3019  (Blu,  B3u)e;  vvs 

CH2  asym.  str. 

2875  (forb.) 

2866  (B3g) 

CH2  asym.  str. 

2874  (x;44) 

2877  (B2u);  vvs 

CHi  sym.  str. 

2868  (z;8) 

2884  (Ag)f 

CHi  sym.  str. 

2868  (y;8) 

2825  (Blu);  vvs 

C=C  sym.  str.g 

1694  (z;7) 

1680  (Ag) 

C=C  asym.  str.8 

1561  (z;5) 

1639  (B3u);  vvs 

CH2  scissor 

1435  (y;15) 

1430  (Blu);  vvs 

CH2  scissor 

HC=CH  i-p  bendd 

1430  (z;2) 

1405  (Blu);  vs 

1428  (Ag) 

HC=CH  i-p  bendd 

1377  0;1) 

1378  (B2g) 

CH2  wag 

1326  (z;3) 

1358  (B3u);  m 

CH2  wag 

1336  (y; 6) 

1280  (B2g)h 

CH2  twist 

1202  (forb.) 

1249  (Au)h 

CH2  twist 

1197  (x;  1 ) 

1242  (Blg)h 

HC=CH  i-p  bendd 
HC=CH  i-p  bendd 

1164  (z;0) 

1159  (B3u);  m 

1200  (Ag) 

"ring  stretch”8 

1050  O’)  19) 

1038  (B2g) 

"ring  stretch”8 

980  (y;  1 1 ) 

887  (B,„) 

Calculated  frequencies  are  (B3LYP/6-31G*:PM3)  results,  scaled  by  a  factor  of  0.9613  (Ref.  [41]).  Only  calculated  modes  above  ~1000 
cm-1  are  given.  For  IR  and  Raman  data,  symmetry  assignments  for  the  D2h  point  group  are  given. 

a  str.  =  stretch;  sym.  =  symmetric  (i.e.,  in-phase);  asym.  =  antisymmetric  (i.e.,  out-of-phase);  i-p  =  in-plane.  The  mode  descriptions 
(“sym.  str.",  etc.)  are  those  given  in  [37,40]  for  the  free  molecule. 

b Calculated  (this  work),  “x”,  “y”  and  “z”  indicate  polarization  behavior  in  IR  absorption  (axes  defined  in  Fig.  4b).  “Forb.”  in¬ 
dicates  "IR-forbidden”  for  C2v  symmetry.  The  number  following  the  polarization  label  is  the  computed  relative  intensity  on  a  scale  of 
0-100. 

0  Relative  intensities  [39]  are  indicated  qualitatively:  vvs = very  very  strong;  vs  =  very  strong;  m  =  moderately  strong. 
d  There  is  only  one  FI— C=C— H  group  in  the  adsorbate;  hence,  only  one  mode  remains  from  each  pair  of  modes  in  the  free  molecule. 
“These  modes  are  accidentally  degenerate  in  the  free  molecule. 

rThis  mode  is  believed  to  be  shifted  to  higher  energy  by  a  Fermi  resonance  [40]  with  the  overtone  of  the  CH2  scissor  mode. 

8 In  the  adsorbate,  the  two  C=C  bonds  are  no  longer  equivalent.  The  higher-  (lower-)  energy  C=C  stretching  mode  involves  mainly 
the  H— C=C— H  (Si — C=C— Si)  group.  For  the  same  reason,  the  “ring-stretching”  modes  differ  from  those  of  the  free  molecule. 

h  These  are  estimated  frequencies  [40],  The  A„  modes  are  IR-  and  Raman-forbidden  in  the  D2h  point  group.  The  other  estimates  are 
for  modes  which,  although  allowed,  are  too  weak  to  be  observed. 

planes  and  z  as  the  twofold  rotation  axis.  3  Upon 
reduction  from  D2h  to  C2v,  all  modes  become  IR- 
allowed  except  those  belonging  to  the  B3g  and  Au 
representations  of  D2h.  These  correlate  [43]  with 
the  A2  representation  of  C2v  which  is  IR-forbid¬ 
den.  The  Ag  and  B3u  modes  of  the  free  molecule 
(which  correlate  with  A!  of  C2v)  are  allowed  in  z- 
polarization.  The  Big  and  B2u  modes  (which  cor¬ 
relate  with  Bi  of  C2v)  are  allowed  in  x-polarization. 


3  In  the  convention  used  here,  the  x'-axis  of  D?h  is  preserved 
as  the  twofold  rotation  axis  upon  reduction  to  C2y 


The  B2g  and  Bhl  modes  (which  correlate  with  B2  of 
C2v)  are  allowed  in  ^-polarization.  If  the  experi¬ 
ment  is  done  on  a  two-domain  surface,  with  equal 
numbers  of  rings  lying  in  the  (pc,  z)  and  ( v,  z)  planes, 
no  distinction  can  be  made  between  x-  and  y- 
polarized  modes.  Of  course,  a  polarized  MIR 
experiment  on  a  single-domain  sample  would 
provide  more  detailed  data.  Such  a  single-domain 
c(2x2)  structure  has  been  used  previously  in  sur¬ 
face-sensitive,  polarized  x-ray  absorption  experi¬ 
ments  [44]. 

In  Table  4,  the  polarization  of  each  calculated 
mode  was  determined  by  examining  the  atomic 
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displacements  and  applying  the  C2v  symmetry 
operations.  The  polarization  assignments  were 
thus  determined  independently  of  the  D2h  <->  C2v 
correlations  given  above.  The  results  show  that  of 
all  the  adsorbate  modes  above  1000  cm-1  only  two 
are  x-polarized.  In  addition  to  the  adsorbate  vi¬ 
brations,  the  calculation  gave  a  — C=C—  stretch¬ 
ing  mode  at  2028  cm-1  on  the  bare  surface  which  is 
allowed  in  z- polarization,  given  the  C2v  symmetry 
of  the  c(2x2)  surface  unit  cell  (Fig.  1).  The  ab¬ 
sorption  is  predicted  to  be  relatively  strong  (48  on 
the  scale  of  Table  4)  and  is  comparable  in  fre¬ 
quency  to  the  corresponding  mode  (2132  cnr1 
[45])  in  the  Raman  spectrum  of  liquid  H3SiC= 
CSiH3 .  The  removal  of  this  mode  upon  adsorption 
on  the  c(2x2)  surface  may  be  detectable  in  an 
MIR  experiment,  and  observation  of  such  a  high- 
frequency  stretching  vibration  would  be  further 
proof  of  the  triple-bond  nature  of  the  — C=C— 
bridge. 

2.5.  The  transition  state 

Fig.  5  shows  an  energy-level  diagram,  following 
Okamoto  [11],  which  compares  the  [2  +  4]  cyclo¬ 
addition  reaction  of  1,3-butadiene  with  molecular 
C2H2  and  with  the  — C=C—  bridge  on  the  c(2x2) 
surface.  The  energies  of  the  molecular  HOMOs 
were  obtained  from  photoemission  data  [46]  and 
the  HOMO-LUMO  separations  from  B3LYP/6- 
31G*  calculations  for  the  respective  free  molecules. 
The  experimental  [46]  1,3-butadiene  HOMO  en- 


Fig.  5.  Schematic  energy  level  diagram  showing  the  alignment 
of  the  1,3-butadiene  HOMO  and  LUMO  with  those  of  free 
acetylene  (C2H2)  and  with  the  — C=C—  n  and  ir*  orbitals  on  the 
SiC(l  0  0)-c(2x2)  surface. 


ergy  is  that  of  the  stable  trans  form.  The  position  of 
the  SiC  conduction  band  minimum  (CBM)  below 
vacuum  is  given  by  the  measured  [47]  electron  af¬ 
finity  for  the  c(2x2)  surface,  4.4  eV.  The  valence 
band  maximum  (VBM)  lies  lower  than  the  CBM  by 
the  cubic  SiC  band  gap  energy  [48]  of  2.40  eV.  The 
positions  of  the  — C=C—  n  and  n*  states  relative  to 
the  VBM  and  CBM  were  obtained  from  the  theo¬ 
retical  results  of  Sabisch  et  al.  [27]. 

For  C2H2,  the  large  energy  difference  between 
the  LUMO  of  the  “dienophile”  (C2H2)  and  the 
HOMO  of  the  diene  means  that  the  reaction  will 
be  slow  under  normal  conditions.  In  fact,  it  has 
been  observed  only  in  the  presence  of  metal  cata¬ 
lysts  [49,50].  However,  the  dienophile  HOMO  on 
the  c(2  x  2)  surface  is  at  a  much  higher  energy  than 
in  free  C2H2  which  reduces  the  energy  difference 
for  the  so-called  “reverse-electron-demand”  reac¬ 
tion,  in  analogy  with  the  C(100)-(2xl)  surface 
[11],  The  energy  difference  of  5.18  eV  (Fig.  5)  is 
somewhat  greater  than  the  corresponding  value  of 
4.60  eV  computed  [11]  for  this  reaction  on  the 
C(1  00)-(2x  1)  surface  modeled  as  a  C9H12  cluster. 

Fig.  6  shows  the  transition  state  obtained  in  a 
QST3  optimization  [21].  Because  of  the  computa¬ 
tional  effort  involved,  only  (B3LYP/6-31G*:UFF) 
calculations  were  done.  A  UB3LYP  calculation 
showed  no  unpaired  spin  density  and  gave  the 
same  optimized  energy.  As  noted  above,  the  results 
are  expected  to  be  essentially  the  same  as  those 
obtained  with  RHF/PM3  as  the  low-level  method. 
The  appearance  of  a  single  imaginary  frequency 
shows  that  the  structure  corresponds  to  a  first- 
order  saddle  point.  Bond  lengths,  defined  as  in  Fig. 
4b,  are  given  in  Table  2.  Of  most  interest  are 
d4  =  2.43  A  (the  length  of  the  bond  forming  be¬ 
tween  surface  and  adsorbate  C  atoms)  and  the 
angle  of  cj)  =  51.8°  between  the  1,3-butadiene 
plane  and  the  surface  normal  (cf.  Fig.  6b).  These 
are  similar  to  the  corresponding  values  of  2.64- 
2.68  A  [4,11]  and  49°  [4]  for  C(1  00)-(2x  1). 

The  activation  energy,  A EA  =  +6.1  kcal/mol,  is 
somewhat  greater  than  for  the  same  reaction  on 
C(100)-(2xl)  [4,11].  It  is  useful  at  this  point  to 
consider  whether  it  is  reasonable  to  expect  the  re¬ 
action  studied  here  to  occur  under  conditions  ac¬ 
cessible  in  a  typical  ultra-high  vacuum  (UHV) 
chamber.  Experimentally,  the  C(  1  0  0)-(2  x  1 )  re- 
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(b) 


Fig.  6.  Schematic  diagram  showing  part  of  the  duster,  together 
with  1,3-butadiene,  in  the  transition-state  configuration,  (a)  and 
(b)  are  viewed  along  the  x-  andy-axes  respectively  (cf.  Fig.  4b). 


3.  Summary  and  conclusions 

The  reaction  of  1, 3-butadiene  with  SiC(lOO)- 
c(2  x  2)  surface  has  been  considered  using  compu¬ 
tational  methods.  The  results  are  as  follows: 

(1)  Semiempirical  and  MM  calculations  predict 
that  the  [2  +  4]  cycloaddition  product  with 
the  — C=C—  bridge  is  the  most  stable  of  six 
hypothetical  adsorbate  structures. 

(2)  A  cluster  model  for  the  bare  c(2x2)  surface 
was  constructed  using  ME  with  divalent  link 
atoms  [24].  A  hybrid  calculation  with 
B3LYP/6-31G*  as  the  high-level  method  and 
either  semiempirical  (RHF/PM3)  or  MM 
(UFF)  methods  in  the  low-level  gave  an  opti¬ 
mized  structure  that  compares  well  with  previ¬ 
ous,  fully  ab  initio  results. 

(3)  The  cycloaddition  product  was  treated  by  in¬ 
cluding  the  1,3-butadiene  adsorbate  in  the 
high-level  section.  For  the  [2  +  4]  reaction,  the 
resulting  structure  resembles  1,4-cyclohexadi- 
ene.  A  reaction  energy  of  A ER  =  —50.3  kcal/ 
mol,  with  an  activation  energy  of  A EA  = 
+6.1  kcal/mol,  was  obtained.  The  correspond¬ 
ing  [2  +  2]  reaction  product  was  found  to  lie 
23.5  kcal/mol  higher  in  energy. 

(4)  The  vibrational  spectrum  of  the  adsorbate  was 
analyzed.  Normal-mode  frequencies,  absorp¬ 
tion  intensities  and  polarization  effects  ex¬ 
pected  in  IR  spectroscopy  were  described. 


action  is  found  [51-53]  to  be  facile  at  room  tem¬ 
perature,  and  chemisorption  saturates  at  a  1,3- 
butadiene  exposure  of  ~30  L  (1  Langmuir 
(L)  =  1(V6  Torrs)  [54].  Assuming  similar  pre¬ 
exponential  factors  for  the  two  systems,  A EA  w  3 
kcal/mol  for  C(100)-(2xl)  and  ~6  kcal/mol  for 
SiC(l  0  0)-c(2x2),  as  in  Table  3,  suggests  that  the 
reaction  will  be  slower  by  a  factor  of  ~  1 50  for  the 
c(2x2)  at  room  temperature.  The  required  expo¬ 
sure  of  ~5  x  103  L  for  SiC(  1  0  0)-c(2  x  2)  is  certainly 
feasible  in  a  UHV  system  for  a  simple  hydrocar¬ 
bon  like  1,3-butadiene.  This  is  especially  true  in  an 
IR  spectroscopy  or  scanning  tunneling  microscopy 
experiment,  for  which  it  is  unnecessary  to  re¬ 
evacuate  all  the  way  back  to  UHV  conditions  after 
exposure  in  order  to  record  data. 
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